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The notion of f-oscillators generalizing q-oscillators is introduced. For classical and 
quantum cases, an interpretation of the f-oscillator is provided as corresponding to a special 
nonlinearity of vibration for which the frequency of oscillation depends on the energy. The 
f-coherent states (nonlinear coherent states) generalizing q-coherent states are constructed. 
Applied to quantum optics, photon distribution function, photon number means, and disper- 
sions are calculated for the f-coherent states as well as the Wigner function and Q-function. 
As an example, it is shown how this nonlinearity may affect the Planck distribution formula. 
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1 Introduction 

In quantum physics, harmonic oscillators are synonymous with creation and annihilation oper- 
ators. For this reason, in the first attempts to realize Hopf algebras (quantum groups) in terms of 
creation and annihilation operators (a generalization of the Jordan-Schwinger map) the resulting 
oscillators were named q-oscillators. This pervading property of the oscillator formalism in many 
physical situations has induced a lot of interest in looking for physical consequences, where honest 
oscillators are replaced by q-deformed ones (partition functions, field theories, nonlinear optics, 
etc.). 

Coherent states, defined through creation and annihilation operators, provide us with a beau- 
tiful connection between quantum and classical oscillators. The notion of coherent states [1-3] 
permitted the use of language and intuition developed from the study of the classical mechanics 
of harmonic oscillators in order to treat their quantum counterpart, because the trajectory of 
the center of quantum coherent packet is the same as the classical trajectory and the width of 
the packet is the minimal possible one in the frame of Heisenberg uncertainty relation Q. The 
notion of coherent state turned out to be appropriate also to describe simple quantum systems like 



spin and cyclotron motion of a charge in magnetic field ||. On the other hand, the notion of 
the quantum q-oscillator |7|, || was interpreted [|, [K| as a nonlinear oscillator with a very specific 



type of nonlinearity, in which the frequency of vibration depends on the energy of these vibrations 
through the hyperbolic cosine function containing a parameter of nonlinearity. 

This interpretation of q-oscillators becomes obvious if one used the classical counterpart of the 
original quantum q-oscillators. This observation suggests that there might exist other types of 
nonlinearity for which the frequency of oscillation varies with the amplitude in a manner different 
from the cos/i-dependence; we will label this dependence by a function /. Such classical oscillators 
(and their quantum partners) may be called f-oscillators |J. It is interesting to consider the 
statistical mechanics of a gas of deformed oscillators (free energy, partition function behaviour) 
and to compare it with the one associated with standard oscillators. 

The problems have obvious counterparts in quantum mechanics. Here, the role of the phase 
diagram is played by the eigenstates of the Hamiltonian. For stationary systems, one could consider 
such changes of the Hamiltonian, which is an integral of motion, that produce new Hamiltonian 
which is some function of the initial one. Then if there are no degeneracies in the spectra for 
the initial Hamiltonian, the eigenstates of the new Hamiltonian coincide with the old ones. But 
for the new system the energy spectra are different. This produces time evolution of the phase 
factors of the eigenstates such that these vary with different velocities in complete analogy to the 
classical motion of the corresponding deformed classical systems, moving along their trajectories 
in phase space with reparametrized velocities. 

In fact, in more general situations the new quantum systems having the same stationary eigen- 
states as the initial ones possess the new Hamiltonian which is a function of the usually commuting 
time independent integrals of motion (complete set of observables) . This is the analogue of the 
new classical Hamiltonian which was deformed using reparametrizations depending on all available 
classical invariants which are in involution. 

The aim of our paper is to study the behaviour of classical and quantum systems belonging 
to the subclass described above, and to clarify the role of nonlinearities corresponding to these 
systems. This goal is motivated by the fact that q-oscillator belongs to the system of the subclass 
corresponding to the specific q-nonlinearity ||. 

In classical case, we consider simple linear systems (oscillators) and their deformations pro- 
ducing nonlinear integrable systems as well as symmetries based on the nonlinear noncanonical 
transform of the conjugate variables preserving the vectorfield. In the quantum counterpart, we 
study such systems which differ in Hamiltonian but have the same set of eigenstates. In these 
cases, we analyze the possibility of extending the notion of coherent states of usual harmonic 
oscillator to the case of f-oscillators. Algebraic extensions of the notion of q-oscillator coherent 



states have been discussed in [11, 12] and applications in [13]. The particular case of f-coherent 



states called also as nonlinear coherent states for the function / expressed in terms of Laguerre 
polynomials was shown to be created for trapped ion in |fH| . Shortly f-oscillators were discussed 
in p). 

We will study some physical consequences of the existence of f-coherent states like the change 
of the particle distribution function, the possibility of having super- or sub-Poissonian statistics, 
influence of f-nonlinearity in the black body radiation formula with the particular example of 
the q-oscillators. As a particular example we will apply these results to q-nonlinear systems (q- 
oscillators) and will show that for q-coherent states there exists sub-Poissonian statistics which 



mean that q-nonlinearity of their fields decreases the fluctuation of the particle number in a 
q-coherent state. 

The next three sections (2, 3 and 4) illustrate in detail the situation in classical mechanics 
with the examples of one-dimensional, two-dimensional and three-dimensional oscillators. The 
introduction of the one- dimensional quantum analogues (the quantum f-oscillator) is given in 
Section 5 and in the following one (Section 6) some algebraic relations are shown for the operators 
describing them. Then the eigenstates of a one-mode f-annihilation operator (f-coherent states) 
are considered in the Fock space (Section 7) as well as in different representations (Wigner and 
Husimi) in Section 8. Properties of such states are studied in Sections 9, 10 and 11, namely, their 
evolution and completeness relations with a remark on the Stone-von Neumann theorem, with the 



main end of underlining that they are not always coherent in the ordinary sense | lq| . Extensions 
to many modes are shown to be possible not only when they are independent of each other but 
also when there is a nonlinear coupling among them (Section 12). The current setup of quantum 
mechanics having been maintained, an application of the above objects in quantum optics is given 
in the last Sections 13, 14 and 15 obtaining new photon distributions, squeezing in the nonlinear 
coherent states, correlation of quadratures, and Planck formula deformation. 



2 Construction of Nonlinear Systems from Linear Systems 

It was shown, that in the classical limit the one-dimensional q-oscillator is represented by a 
reparametrized oscillator ||. The reparametrization is provided by a constant of motion and the 
associated differential equations exhibit a special kind of nonlinearity. Here, we would like to 
consider these systems from a general viewpoint: nonlinear systems as "reparametrized" linear 
ones. 

A linear dynamical system, say on R n , with coordinates x = (xi, x%, . . . , x n ) is described 
for any nxn matrix A by the differential equation 

x j = A k j x k ; AjeR, (1) 

with solutions 

x (t) = e iA x ; x = x (t ) . 

To obtain new (nonlinear) systems from the above, we replace the globally constant matrix A 
with a matrix valued function B (x) and write the equation 

x j = B^)x k (2) 

with 

It means that the matrix elements of the matrix B are integrals of motion for Eq. (^|). For this 
new system, we can write the solutions as 

x(t) = exp [tB (x )] x . 



For each initial condition, Eq. (|2|) reduces to Eq. (|1|) and a particular case is the orbit-dependent 
time reparametrization. 



An example of such a deformed equation in field theory has been shown in [ I7| where a physical 
parameter was made dependent on Cauchy data. A particular family of such systems is when 



5j(x) = /(x)A*, 

with / any constant of motion of the system. For instance, let us consider the three-dimensional 
isotropic harmonic oscillators 

x = cuy; 

(3) 

y = — uy. ; u G R . 

All constants of motion of this system are functions of bij(x % yi — x^y 1 ) and aij{y % yi + x l x^) . 
If we set new frequency 0, = O ( {pij(x l yi — x^y 1 ), aij(xiXj + y%y^) ) and make the replacement 
u — > Q in Eq. (|3]), we get a nonlinear system. Solution to this system is given by 

x i( t ) \ _ ( costtt(x ,y ) smtQ(x ,y ) \ ( Xj(0) \ ,^ 
yj(t) J \ -smtn(x Q ,y ) costn(x ,y ) J \ %(0) J ' 

For each choice of the functional dependence of Q on the constants of motion, we get a foliation 
of the carrier space R 6 with leaves given by 

S A = {(x,y)GR 6 : Q (x,y) = A} . (5) 

For each initial condition x , y G S A , the oscillatory motion has the same frequency for each 
one of the coordinates. 

It is however possible to start with 
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(6) 



and "nonlinearize" the system by making different choices of the constants of motion for the 
different frequencies. For instance, we could take 

= h(xl + yf,xl + yi,xl + yl); 
^2 = f 2 (xl + yl,xl + yl,xl + yl); (7) 
^3 = h{x{ + vl, xl + ylxj + yj). 

Now each mode will be "dynamically coupled" to the others in a different way. 

If T is the dynamical vectorfield, it should be noticed that when we replaced T with fT , i.e., 
we reparametrize T by a constant of motion, we get the same phase portrait for both vectorfields. 
If we reparametrize different modes of vibration differently, we change the phase portrait of the 
new system with respect to the original one of T. Thus, our "nonlinearization" procedure is more 
than a parametrization of the original linear system. 



3 Deformation of Linear Hamiltonian Systems and Sym- 
metries 



(8) 



The above variety of choices has to be restricted if we start with a system admitting a Hamil- 
tonian description and want to preserve the Hamiltonian character after we "reparametrize" it. 
The simplest choice to get a reparametrized Hamiltonian system is the following: Start with a 
linear Hamiltonian system 

/ ±i \ ( dH/d Vi \ 
{Vi ) \ -dH/dx t J ' 

where 

H = XiA*y j + B ij y i y j + C ij x iXj 

(the coordinates yi are the momentum components) and consider the new system associated 
with the Hamiltonian H — f (H) . We get the new nonlinear equations of motion given by 

dyi 

(9) 

OXi 

with f'(H) = df /dH . This system can be explicitly integrated as any linear system can be. This 
has already been demonstrated in |§. 

To make contact with ||, we also show, however, a different path to nonlinear Hamiltonian 
systems obtained from linear systems, from linear symmetries to nonlinear ones. If A is any 
nxn matrix, we denote by Xa the associated vectorfield setting Xa = XiA l -d/dxj . For a 
dynamical system, we prefer using the notation V instead of Xa ■ The Lie algebra of symmetries 
for T contains all linear vectorfields Xb such that [B, A] = ; this follows trivially from 
[Xa, Xb] = X\a,e\ ■ in the case of the isotropic harmonic oscillator, we find that for the Tri- 
dimensional oscillator the symmetry algebra is gl (m, C) . Now we "reparametrize" any one of 
X' B by using constants of motion for T . Had we started with matrix B generating symmetry 
transformations instead of infinitesimal ones, i.e., B~ 1 AB = A (diffeomorphism in the differential 
geometric language), we could generate nonlinear changes of coordinates by following the same 
idea of "reparametrization." Generally, this procedure will turn canonical transformations into 
nonlinear noncanonical transformations. 

We deal with the two-dimensional isotropic harmonic oscillator on R 4 

Xi = Vi] iji = -Xi\ £ = 1,2 (10) 
and define the change of coordinates 

qi = fi(Hi)Xi] 

(11) 

Pi = fi(Hi)yi 



with Hi — x\ + y\ ; : R — > R and no summation on i . In these new coordinates, the 
equations of motion have still a linear form given by 



Pi 



(12) 



Pi 



-Qi 



Without loss of generality, we consider only one degree of freedom. In both coordinate systems, 
the dynamics admits a Hamiltonian description given respectively by 



and 



a) {x,y} = l 



b) {p, q} = 1 ; 



V)\ j t f = i H >f} ( 13 ) 

H = \ (P 2 + q 2 ) ; jj = {H, /} • (14) 

However, using the Poisson bracket of a) to compute the Poisson bracket in b) for (p, q) as 
functions of (x, y) given by the system of equations (0), we find 

{f(H)x,f(H)y} = ^(Hf(H))^l, 

i.e., the nonlinear transformation we have performed is noncanonical. The noncanonical property 
is there even if F(H) is a constant / . To obtain the same right hand side of the Poisson bracket 
in b) , it is necessary to use for the independent variables (x, y) the new Poisson bracket 



{x,y}' 

with it and the Hamiltonian function 

H'(x,y) = 



_d_ 

Jh 



(Hf(H) 



u 2 • ,r) /V + ir) 



(15) 



(16) 



we have another Hamiltonian description in the coordinate (x, y) for V . 
This can be written in symplectic terms as 



and 



Iyoj = —dH ; uj = dx A dy 
d 



(17) 



-dH' ; J 



dH 



(Hf 2 (H))dx Ady. 



Because both H and H are Hamiltonian functions for V , we can use H expressed as function 
of (p, q) through the inverse of ([ll]) and the bracket {p, q} = 1 to introduce a new vectorfield. 
Otherwise, we can use 

-l 



{x,y}' 



d 
dH 



(Hf(H) ) 



and the Hamiltonian 

H=\{x 2 + y 2 ). 

In both cases, we get a "reparametrized" version of harmonic oscillator. In comparing the two 
procedures, it is clear that in one approach we have used the same Poisson Brackets ( "commuta- 
tion relations") as before with a new Hamiltonian which is a function of the standard quadratic 
Hamiltonian. In the second approach, we use the same Hamiltonian function but the "deformed" 
Poisson brackets ("deformed commutation relations"). It is this second viewpoint which allows the 
connection with q-oscillators, while the first one has suggested the interpretation of the classical 
counterpart as a nonlinear f-oscillator. However, it should be clear that both views are legitimate 
and both of them give rise to a nonlinear dynamics out of a linear one by taking recourse to a 
"reparametrization." 

Having in mind the passage to quantum mechanics, we prefer turning now to complex variables 
in R 4 and defining complex coordinates 

a k = 2- 1/2 (x k + iy k ); 

(19) 

a* k = 2- l l 2 {x k -iy k )- k = 1,2. 
The complex structure in R 4 is given by the matrix 



J = 



( 





1 





\ 




-1 






















1 


\ 








-1 


o J 



satisfying J 2 = — 1 . This matrix defines a complex structure commuting with the dynamical 
evolution associated with the isotropic harmonic oscillator 

r = 1 " al ik) 

The equations of motion are 

a = —ia ; 
a = ia . 

For the (p, q) coordinates, we have new complex coordinates 

a = 2- 1/2 {q k + ip k ); 

(22) 

C k = 2- l ' 2 (q k -ip k ) 

and, setting 

n k = a k a* k , 



(20) 



(21) 



we find 



& = fk(n k )a k ; 

(23) 

Ck = fk(n k )a* k . 

This transformation is not "analytic" (we notice that analyticity depends on the complex structure 
and here we have two alternative complex structures compatible with the dynamics). 

The equations of motion for these new variables are the same as ( pi]) as it was seen earlier, i.e., 
Eqs. ([H]) are form invariant under the transformation (Pg). in these complex coordinates, if we 
take the new point of view, stemming from quantum mechanics, which takes the Hamiltonian as a 
primitive concept for the dynamics, we are naturally led to consider the following two Hamiltonians 

Hi = \Y. a ^l (24) 

k 

in the a-coordinates, and 

#2= J £M (25) 
z k 

in the ^-coordinates. To compare, we express them in the same variables to find 

ffi = £$>* ( 26 ) 
A k 

and 

H 2 = WfMn k . (27) 
/ k 

We now use the same bracket for them, say 

{a k , a*} = -iS k j ■ (28) 

We obtain two different dynamical systems, where the one associated with H2 is not even nec- 
essarily isotropic. The evolution goes from a periodic orbit to an orbit whose closure is a two- 
dimensional torus, i.e., the associated systems are completely different. Of course, there is no 
contradiction. Indeed, to have the same dynamics using H2 we should use different Poisson 
brackets, as it has been seen earlier. 

For the nonlinear oscillator obtained by means of the deformation function /, the equations 
of motion are 



a = —i—(nf 2 (n))a: 



a* = i -^-(nf 2 (n)) a* 



(29) 



and are not invariant under the transformation (p3]), which indeed gives the equations of motion 
for another system of our class but with a different deformation function. They admit different 
Hamiltonian descriptions. For instance, 

1. H = aa* ; {a, a*} = i — (nf 2 (n) ) ; lu = ( — (nf 2 (n) ) J da A da*, 

(30) 

2. H = nf 2 (n); {a,a*} = — i; u = daAa*. 

It is now clear, that in the quantum picture these complex coordinates will be replaced by creation 
and annihilation operators. Therefore for the corresponding commutators, we can repeat what we 
have said for the Poisson bracket. 



4 Some Examples 

The first example we consider concerns the classical one- dimensional q-oscillator |1, where 



' sinh Xaa* 
aa* sinh A 



-a 



(31) 



r 



' sinh Xaa* 
aa* sinh A 



a 



The Poisson bracket for the new variables can be computed using (p8| ) and expressed in terms of 
themselves, by the use of the map (£, £*) — > (a, a*) , the inverse of (|23| ) (which in this case exists). 
One then obtains 

{£, H = -2— ^ V /l + l^l 4 ( sinh A ) 2 ' (32) 
so that we can consider a new system described by such variables with Hamiltonian function 



The equations of motion then are 

i = 



A 



sinh A 

(and its complex conjugate) with solutions 

-itX 



l + |£| 4 (sinh A) 2 e 



e(t)=e(0)exp 



sinh A 



1+ | f (0) | 4 (sinh A) 2 



(33) 



(34) 



(35) 



(and its conjugate). Such a system can be rewritten, of course, in terms of (a, a*) variables: in 
these coordinates, the original Poisson bracket is unchanged while the Hamiltonian function is 



(36) 



It is clear, that this new dynamical system has a phase portrait which is the same as the usual 
linear harmonic oscillator. The only difference is in the frequencies, the new one being 



to = — - — - cosh Acta* . (37) 
smh A 

We notice that aa* is a constant of motion for both systems. The deformation function used 
has given both energy and frequency exponentially growing with aa* . To have physically more 
acceptable functions, it is not difficult to consider a slightly different deformation function. 

The other example we consider is the deformation which leads to the classical version of 
harmonious states to be discussed in next sections. Here the deformation function is taken to be 

( 1 \ 1/2 

f( n ) = 1 r(w + i) J ' n = aa *> ( 38 ) 



so that 



The Poisson bracket 



1 





t = tv ; ^ «*• (39) 



where the right hand side is the function of ££* which is its transform by the inverse of (|23f) , 
remarking that the T function on the positive real line has well defined derivative as well as 
inverse. 



5 Nonlinear Oscillators in Quantum Mechanics through 
Noncanonical Transformations 

We wish now to deal with the quantum analogue systems: we will remain in the usual formula- 
tion where the Fock space is better to describe the quantum harmonic oscillator. To go along the 
same lines as in the previous sections, we stay in the Heisenberg picture and write the equations 
of motion for the harmonic oscillator amplitude a 

a = —iioa 

and for its conjugate 

The transformation ( PB| ) is written here as A = a f((va) ; A' = f{a)a) a) . It is noncanonical since 
it does not preserve the commutation relations. The operators A and A^ evolve with the same 
equations, i.e., 

A = -iuA , A ] = iuoA ] 

in complete analogy with the classsical case. We have in our Hilbert space the vacuum state | 0) 
which satisfies 

a | 0) = , 



as well as 

A | 0) = 0. 

This allows us to construct two bases in the vector space having this vector in common. One is 
the standard (Fock) basis 

|n> = ^=|0>, 



which is orthonormal in standard scalar product 

(n | m) = 5 nm . 
Another basis is constructed using the operator A\ 

I n) = -= | 0). 



We define a new scalar product in the same vector space which gives 

(n | m) = 5 nm . 

Provided /(a* a) to be nonsigular, we can then speak of two Hilbert space structures carried by 
the same vector space. 

The adjoint with respect to this new scalar product does not coincide with the old one. We 
can then define the operators 



n) 



\/n + 1 | n + 1) ; b \ n) — ^/n \ n — 1) ; 



where * means the adjoint in the new scalar product. These operators satisfy the commutation 
relations [b, b*] = 1 . Taking the Hamiltonian H = uob*b we have for the operators b, b* the 
equation of motion of the harmonic oscillator. This is actually the situation with A and A^. 

In the second Hilbert space, the operators A and A^ have an identical representation as a 
and a) have in the Fock space and so they satisfy the commutation relation [A, A'] = 1 . Thus, for 
one and the same vector space, we have the possibility to introduce two Hilbert space structures. 
As for the dynamics of the harmonic oscillator, we have two different descriptions, which parallel 
the alternative Hamiltonian descriptions of the classical oscillator. Therefore, much as we did for 
the classical case, we can use the new Hamiltonian and the old commutation relations to get a 
"deformed" dynamics, or vice versa. To keep the current physical interpretation of the operators 
a and a\ we choose to maintain the value of their commutator which is directly connected with 
measurements, while we will consider deformed Hamiltonian operators. Then operator of "energy" 
H is 

H=-(A*A + AAi). (41) 
2 

In the Fock space, its eigenvalues are 



E n = -[(n + l) f(n + l)f*(n + l)+n f(n) f*(n) } 



(42) 



We illustrate the situation with an example and consider the nonlinear noncanonical transforma- 
tion such that 



A = af q (a^a) = a 



I sinh An 
nsinh A 



which is invertible, 



a = A 



In 



N sinh A + V^ 2 sinh 2 A + 1 



where 



XN 



N = A*A. 



1/2 



(43) 



(44) 



(45) 



The operators A, A* acting in the same Hilbert space as the operators a, (the original Fock 
space) satisfy in this case the commutation relations 



[A, A f ] = iV(cosh A - 1) + \/iV 2 sinh 2 A + 1 , 



(46) 



as can be seen expressing all the operators in terms of matrices. 

If one has the linear dynamics for the operator a with frequency equal to unity, i.e., 

a + ia = , 



(47) 



and boson commutation relation for the operators a and a\ the same dynamics exists for the 
operator A : 

A + %A = , (48) 

since 

A + iA = (a + ia)f q (a^a) (49) 

and the function f q is an integral of motion. The Hamiltonian for this dynamics may be taken 
as 

I - . I 

(50) 



H= - In 
A 



AM sinh A + J{AU) 2 sinh 2 A + 1 



1 

2' 



and one obtains 

[A,H] = A. (51) 

For the same dynamics, it is possible however to have a different Hamiltonian formulation. In 
one case, it is related to the above Hamiltonian, while in another Hilbert space we define the 
Hamiltonian 

(52) 



H' = A^A + - 
2 



with commuation relation 
and again 



[A, ^] = 1 , 
[A, H']=A. 



(53) 
(54) 



Thus, we see that analogous to classical mechanics there are possible alternative descriptions of 
a quantum system. For the same equations of motion, for the operators we have two different 
Hamiltonians with corresponding different commutation relations. 

The same situation is anyhow present also for the class of nonlinear oscillators we have consid- 
ered. There, continuing with the example of the f q deforming function, we start with the dynamics 
for the q-oscillator |9j given by the equation 

a + i [ sinh X(a^a + 2) — sinh Xa)a ] a = . (55) 

2 sinh A 

It worth remarking at this point, that if we multiply the last equation from the right hand side 
by the same function f q (a'a) , for instance, we are led to a new f-oscillator. Since this function 
is an integral of motion for the above q-nonlinear equation, we have in fact for operator (|33|) the 
equation of motion 

• sinh A . , , + , . . . 

A + i — - — (cosh Xa'a)A = . (56) 
A 

But after (||), 



a) a = — In 
A 



A^Asinh X + J (At A) 2 sinh 2 X + 1 



the obtained dynamics is different from the initial one. 

We return now to consider Eq. fl5"5j) and the first Hamiltonian description is given by 



(57) 



H = [ sinh A(a^a + 1) + sinh Xa)a ] ; [a, a*] = 1 . (58) 

2 sinh A 

In another Hilbert space, let the operators B and evolve with the equation 

B + i l —— \ sinh X(B ] B + 2) - sinh XB^B 15 = 0. (59) 

2 sinh A J v ; 

If we take the commutation relation 



[B, B f ] = 5 f 5(cosh A - 1) + J(BW) 2 sinh 2 X + 1 , (60) 



the form of the Hamiltonian for this system differs from the form of Hamiltonian 

The important physical consequence of the existence of the same dynamics for quadratures 
with different commutation relations is the possibility of existing identical harmonic vibrations 
of two kinds. One vibrational process respects the Heisenberg uncertainty relation since quadra- 
tures satisfy standard boson commutation relations. Another vibrational process is compatible 
with different uncertainty relation for its quadratures since they satisfy different commutation 
relations. Nevertheless, from the view point of dynamics (harmonic vibrations) both cases are 
undistinguishable. 



6 f-oscillator Operators 



The operators A and A* represent the dynamical variables to be associated with the 
quantum f-oscillators. The well known q-oscillator operators belong to this class. In this section, 



after discussing some algebraic features useful for the f-generalization, we refer also to other 
examples of f-oscillators already taken into consideration. 

We start by recalling some notions about the harmonic oscillator operators a and a' whose 
algebraic structure is [a, at] — 1 . In the Fock space with a = (at)t ; n = a^a , the basis is given 
by the eigenfunctions of n 

n | n) = n | n) ; n £ Z + . (61) 



We have also 



and V / : Z 



^2\n)(n\; (n\m) = 5 nm (62) 



f(n)=J2f(j)\j)(j\ ■ (63) 

3=0 

Consider now a "distortion" of a and a) of the form 

A = af(h) = /(n+l)a; 

(64) 

^ = /t(^) a t = fl t/t(£ + i) 

and note that 

[A,n]=A; [^,71] = (65) 

The functions we are considering can be made dependent in general, also on continuous parameters, 
in such a way that for particular values of them the usual annihilation and creation operators are 
reconstructed. We will say then that we are in presence of continuous deformations. This was the 
case of q-deformations [18]. In principle, one may consider discontinuous deformations, too. 
Since 

oo 

a = v 7 ^ I n ~~ -0 ( n I ' 

71=0 

(66) 

oo 

a' — \fn \n)(n — l \ , 

n=0 

the same Fock space is a carrier space for A and A\ i.e., 

oo 

A = Z! Vnf(n) | n - l)(n | ; 

n=0 

(67) 

oo 

At = Vnf*(n) | n)(n-l | . 



This realization may or may not be irreducible depending on the assumed functional form of 
/(»)■ 



Following the choice of not deforming the commutators of the physical variables (a, a') , the 
commutator between A and can be easily computed and by using fl67|) reads 



F(n) = [A, At] = (n + 1) f(n + 1) ( /(n + 1) ) f - n /(n) ( /(n) )t , (68) 

while the q-commutator is 

G(n) = [A, A\ = AA^ - qA^A 

= (n+l)/(n + l)(/(n + l)) t -gn/(n)(/(n)) t ; g e (0, 1] . (69) 

We have, also 

F(n) - G{n) = (g - 1) n/(n) / f (n) . (70) 

Since only ff* occurs, the phase of / is irrelevant and we may, without loss of generality, choose 
/ to be real and nonnegative: 

/t(n) = /(n). (71) 

Alternately, given the functions F or G, assumed hermitian, we obtain the following solution to 
Eqs. (0) and (§§ 

n— 1 x 



/( n ) = 4r lE^O')] ^ ( 72 ) 



?? 



and 



'n— 1 x 



/(n) = 4=|E9 , ' G ( n -J'- 1 )| ; n ^°> ( 73 ) 



respectively, with / (0) arbitrary in both cases. Such solutions are unique, having been obtained 
by construction. Of course, we may use complex / to construct the functions F and G in 
( |68| ) and (|69]) . However, to obtain / from F and G, we remark that F has to be real as 
well as G when q is real. Then in (|6~4D , we have the freedom of choosing an n dependent 



phase of /, which corresponds to construction of generalized coherent states [19, g(J. These 
generalized coherent states were analyzed in were dyfferent types of interesting new states 
were introduced. 

As was mentioned earlier, in the case of the q-oscillator operators, the function / depends 
also on a continuous parameter in order to obtain the harmonic oscillator operators as a limiting 
case. Starting with the q-commutation relation [0, [8] 

G{n) = q- fl - A = lng; AgR, (74) 



using (f73[) one obtains 



ft \ h-q n -q~ n / sinh An . 

f( n ) = \ ~ = V n-r = / ff (n), 75 

V n g — q | n smh A 



setting 

/,(0) = 1. (76) 



We see from (|58|) that the eigenvalues of the "energy" operator grow exponentially with increasing 
occupation number. 

One can remark that the phase operators V, [[22] are actually deformations of the Bose 



operators of the kind we are studying and lead to the harmonious states |23| to be considered 
below. In this case, 

/(n) = 4= = /fc, (77) 



so that 

A\n)=\n-1); A f | n) =| n + 1} ; n^O; A | 0) = . (78) 

When many degrees of freedom are involved, we have two possible choices. The first one 
defines q-oscillators which satisfy q-commutation relations among a single degree of freedom, but 
mutually commuting between different degrees of freedom, possibly with different functions for 
the various degrees of freedom. Another choice is to make / dependent on the total occupation 
number operator 

"tot = Yl ™* ■ (79) 

i 

For two degrees of freedom, this was discussed in M. 



7 Nonlinear Coherent States 

Coherent states were originally introduced as eigenstates of the annihilation operator for the 
harmonic oscillator and then widely used in physics, particularly in quantum optics. This is 
therefore a concept of algebraic origin and having now constructed a similar annihilation operator 
it is natural, following the same procedure, to construct a new class of f-coherent states in the 
Fock space. This construction is in general different from other ones |]24j|. Further the f-coherent 



states may not be preserved under time evolution. Nevertheless, we are willing to call them f- 
coherent states for an easy identification, of the kind already proposed for the eigenstates of the 
q-annihilation operator which were named q-coherent states J7|, f25| . 

Let us take for the one-mode case the operator A ( p4[) . Then one can consider the eigen- 
functions | a, f) of A in a Hilbert space. They therefore satisfy the equation 

A | a, f) = a | a, f) ; aeC. (80) 

Looking for the decomposition of | a, f) in the Fock space 



a 



n=0 



we obtain for the coefficients c n the following recurrence relation 



c n +Wn + 1 f(n + 1) = c n a . (82) 



This gives 



in which 

[/(n)]! = /(0)/(l)---/(n)- (84) 

To fix Co , we use the condition 

(a, f\a,f) = l (85) 

and obtain 

/ oo i |2n \-V2 

Co= lS»wj ■ (86) 

To emphasize the dependence of Co on / and \a\ , we will write 

co = N La (87) 

and in order to have states belonging to the Fock space it is required that 

0<iV /ia <oo, (88) 

therefore not any / and |a| are allowed. We will denote with p the positive number such that, 
given /, the above series converge V|a| < p. The scalar product is easily written 



<" I P) = N f , a N f>p ) (a*f3r ! M, \P\ < P- (89) 

n=0 n - I U \ n )\- 



No further constraints are then put on / and p. 

It should be remarked furthermore that, given C (n) = C n any real function on Z + , the 
state | a, C) defined by 

oo 

\a ) C) = Y J C n a n \n) (90) 

n=0 

is an eigenfunction of some A . In fact, the corresponding function / is found to be 

! ^Tn C ~^- <91) 
Such eigenstate can be normalized if the / so obtained satisfies (p6|). In the case 

/(n) = l, (92) 

| a, 1) denotes the usual coherent state and 

iVi,a = expf-tfj , (93) 

a can be any complex number. 

As anticipated, the known q-coherent states [7|, turn out to be a particular case of f- 
coherent states, which we might also call f^-coherent states. Normalization factor of such states 
is 



in which 

sinh An sinh X(n — 1) 
sinh A sinh A 
Using the notation (|84]) we can also write 

[n}\ = [n f q {n)]\ . (96) 

It is seen that a can be any complex number. For the scalar product, we have 

/ 00 I „,l 2ri \ / co \ o\2n\ oo 



(a 



Harmonious states |23] are eigenstates of the annihilation operator deformed by the factor 
fh (|77|) , to which corresponds the normalization 

N h>a = (l-\a\ 2 )- 1/2 , (98) 



and the acceptable a must have modulus less the 1 . Following ( |30"D they can be denoted | a, fh) 
and their scalar product is 

<«, /* 1 & /*> = (j3|^J (r^wj (") 

8 Nonlinear Coherent States in Different Representations 

Since the state | a, f) is given as series of Fock states, we can easily write the wave function 
of these states in different representations explicitly. 
In coordinate representation, the wave function is 

«=^.^Ey imw H - lxh (100) 

where H n is the Hermite polynomial of degree n . 

For the momentum representation, the formula is the same. 

For the Bargmann representation (the usual coherent states), the wave function (z | a, f) , 
where we use the basis | z) [z G C) with a \ z) = z \ z) , takes the form 

^> =N °->^U z ' ar wm- (101) 

For a continuous parameter / , in the limit / — > 1 the usual wave function is recovered 

V£\ = exp f |z|2 ~ |a|2 + J . (102) 



In the Wigner-Moyal representation |26|, the density matrix for the f-coherent state reads 



oo oo 



1 



1 



o^ m\ [f(m)}\[f(n)} 



■a n (-a*) m [V2(x - ip)} m ~ n L™- n {2{x 2 + p 2 )) , 

(103) 



where denotes a generalized Laguerre polynomial. For the particular case of q-oscillator 

f = fq i 



|„, |2n\ 1 oo oo i I 



[n}\ 



E E T^TTT^T"^-"*)" 1 ^ - ^)] m "^r n (2(x 2 + p 2 )) 



m=0 



^ (104) 

Finally, we consider Husimi-Kano |2^J Q-function of f-coherent states. In the definition, 
Qtp(z,z*) is the diagonal matrix element of the density operator | ip)(ip | for the state | ip) 
in the usual coherent state basis. For an f-coherent state | a, f) , we can write 



oo oo ( ~* ~,\m 



Q f (z, ^He-^X.EE 



(z*a) m (za*Y 



m=0 



^ m\ [f(m)]\n\ \f(n)]\' 



(105) 



For q-coherent state, the Q-function is 



oo \ a \2n\ z oo oo 



EE 



(z*a) m (za*) n 



a* W / m =on=o ^/m![m]! ^n![n]! 
For the harmonious states, the Wigner-Moyal function is 



(106) 



W h (x, p) = 2 (1 - |a| 2 )e-^ 2+ P 2 ) E E \/^T ""(-aTfv 7 ^ - zp)] m ""Lr n (2(x 2 + p 2 )) (107) 



m=0 n=0 



and the Husimi-Kano function 



Q h (z, z*) = e-^(l-|a| 2 )EE 



00 00 (2*a) m (^a*) n 



=0n=0 Vm! V™ 



'108) 



9 f-Coherent State Evolution 

Here we offer a few considerations with the aim of underlining the peculiarities of the new 
classes of states in the Fock space, cautioning against relying too much on the abuse of their 
name. 

We consider our field mode evolution to be guided by the equation of motion with the quantum 
Hamiltonian 

H = ~(AA + AAi), (109) 
i.e., in the variables (a, a*) it is the system with the selfinteraction described by the Hamiltonian 

H = \[nf 2 {h) + {n + l)f 2 {n + l)}. (110) 



Then the evolution operator 

U(t)=exp(-itH(n)) (111) 

for this quantum nonlinear oscillator gives the following solution to the Heisenberg equation of 
motion for the operator a (t) 

a(t) = U ] {t)aU{i) = a exp [-iu(n)t] , (112) 

where 

u (n) = \ [{n + 1) f(n + 1) - (n - 1) f\n - 1)] . (113) 

Thus, we see then that also the quantum f-oscillator vibrates with a frequency depending on the 
amplitude. 

Turning to the Schrodinger picture, we can remark that at time t the harmonic oscillator 
has become a deformed oscillator of the kind we discuss. The deformation function is actually 
complex of modulus 1 and we will denote it with 



Ff(n, t) = exp 



-it 



(n + 1) f 2 (n + 1) - (n - 1) f 2 (n - 1) 



(114) 



It is possible, in fact, to introduce the notion of f-coherent states also for complex deformation 
functions / as all formulae go through unaltered. 

Then it can be seen, that if initially the state was the usual coherent state, i.e., in an eigenstate 
of the operator a , then it evolves becoming at time t an -coherent state. Physically it means 
that f-nonlinearity creates the Ff(t) -coherent states in the evolution of a usual coherent state. 
Due to this, the photon statistics of the initial coherent state, to be discussed later, is influenced 
by the f-nonlinearity of the field vibrations. It is evidently different from the usual coherent states. 
Interesting physical example of the f-nonlinear systems is quartic nonlinear oscillator usually used 
for modelling the Kerr medium. 



10 Irreducibility and Deformation 

The usual Stone-von Neumann theorem states that the operators q and p (or a and a) ) 
have no invariant subspaces in the Hilbert space of the oscillator states. If / (n) is chosen to 
have no zeroes in Z + , the operators A and are irreducible over the Fock space. If there 
are one or more double zeroes, the Fock space breaks up into a countable number of irreducible 



representations (compare with Master Analytic representations p8| ). If the zeroes are simple 
zeroes, some of the reduced pieces will not allow a unitary resepresentation. 

It is easy to prove, that if the function / (n) has no zeroes at positive integers, the Stone-von 
Neumann theorem can be extended to the case of the operators A, A* . So for the q-oscillator, 
we are just in this case. Here the map a — > A is invertible and the statement is obvious. More 
interesting situations arise when, for example, the function / (n) has one double zero at the 
integer n , i.e., / (no) = f'(n ) = . Then the subspace of the states 



no 

^> = H s n \n) (115) 

n=0 



is an invariant subspace for the operators A, constructed by means of such function / (n) . 
The subspace of the states 

oo 

|V')= E S n\n) (116) 
n=no+l 

is another invariant subspace for the above operators. Thus, the coherent states defined in this 
case do not contain the states with photon numbers less or equal to n . In this case, the coherent 
state contains the states with photon numbers starting from n + 1 , 

!«,/) = N E ~7=T7TTxu I n >> ( 117 ) 



[/(")]« 

where 

n! = n(n-l)(n-2)---(7io + l), (118) 
[/(n)]! = /(n)/(n-l)..-/(no + l), (119) 

and 

(\ -1/2 
00 Id 2 ™ \ 

n=n +l ™ ! ( [/(")]!) / 



11 Completeness Relations 

In this section, we will show how the f-coherent states form a complete system of states in the 
Hilbert space for nonsinguar / (n) , so that any state vector may be represented as a superposition 
of the f-coherent states. There may be different forms of completeness relations since the set of 
f-coherent states are over complete. 

We introduce first an integral representation for the identity operator which uses the analyticity 
of the f-coherent states and the Cauchy theorem. By construction, the state N7 \ a) is an 
analytic vector valued function of the complex variable a . Hence, the following relation holds 

\n) = ^r[f(n)]\jdaNjX \ a, f) a— 1 . (121) 
Inserting this formula into the known resolution of identity QB2"D we obtain 

1 = --Lf^<f cfdadp* \a,f)((3,f\ NjXNj^Wy n - l n\{[fXn)]\y . (122) 

These are line integrals along contours taken in a region of the complex planes where the con- 
vergence is guaranteed for the series considered in Section 7 to define the normalization constant 
Nf t a . The introduced nondiagonal resolution of identity permits us to calculate the coefficients 
necessary to represent any vector as a superposition of f-coherent states. We also point out that 
the components of the f-coherent state in any basis are the generating functions for components of 
the Fock states. It means that for any matrix representing an operator in Fock basis, the matrix 
elements of the same operator in f-coherent state basis are the generating functions. 



In order that such states can be considered as coherent states in the usual definition |]2"4]| , one 
should write a diagonal resolution of the identity 



dfi(a) | a, f)(a, f | , (123) 
where \x (a) is the weight function. Then the following relations have to be satisfied 

2tt / V+ 1 [N f (p) ]V (p) dp = n\ ( [f(n))\ f ; G Z+ . (124) 



These are actually an infinity of moment equations for the measure \x . For the usual coherent 
states, as well as for the eigenstates of the q-deformed annihilation operators, such measures exist 



25| and in both cases the integral is over the entire complex plane. 

As far as the harmonious states of Section 7 are concerned, applying the general resolution of 
identity (0), one obtains 

da d(3* 1 



4vr 2 



<*p - 1 ^/(i-M 2 )(i-i/?n 



<x,h)(P,h\ (125) 



after having inserted the normalization constants ( |98|) for the harmonious states | a, fh) and 
| (3, fh) ■ We can make this integral along two contours as an integral over the phases of a and 
(3 , if the countours are chosen according to 

a = ae^- f3 = ae^, (126) 

and < a < 1 . Then the double contour integral transforms into 

fl 2 pit r2n £«(<£-' */0 



We conclude remarking that for each a G (0, 1) there is a completeness relation in terms of 
projectors on harmonious states. 

It is possible, however, to have a resolution of the identity with the integral depending on one 
parameter only, once states with norm not necessarily strictly positive are allowed. In p3[ , such 
states have been considered after having introduced the following scalar product 



Then, since 



we have 



(a | p) = (1 - a*P) . 

oo 

a, f h ) = (l-\a\ 2 y 1/2 X> n |n>; |a| < 1 



n=0 

| | a | e if)(| a | e <« | ^ = V |a| 2n (1 - M 2 )- 1/2 | n)(n 

27T q 



and 



Hence, 



JO oo 

\ a \-i e * 9 }(\a\ e ie | ^ = (1 - |«| 2 )- 1/2 (1 - |«r 2 )- 1/2 V | n)(n 



(2 - \ a \ 2 - |a|- 2 )V2 / | | a |-i e je )(|a|e 





d0 _ 
27T ~ 



12 Two-mode f-Coherent States 



For usual multimode harmonic oscillators, there exist generalized correlated states p9|, |30 



in 

which the mode quadratures are statistically dependent. The quasidistributions for these states 
have Gaussian form and the photon distribution function is described by multivariable Hermite 
polynomials . It is possible to extend nontrivially the construction of one- mode f-oscillator to 
many modes. In particular, for two-mode state we consider the two operators 

A i = a i f{n); i = 1, 2, (128) 

where 

h = hi + h 2 ; hi = a\di (129) 

(the operators satisfy boson commutation relation). These operators commute and for this 
reason we can construct algebraically the two- mode f-coherent state | ai, a 2 , f) defined by the 
following equations 

Ai \ an, a 2 , f) = ai \ ai, a 2 , f) , % = 1, 2 . (130) 
Considering the series expansion 



oo oo 



I Oil, ot 2 , f) = E E c «i,«2 \ni,n 2 ), (131) 

m=0 712=0 

where the Fock states | ni, n 2 ) satisfy 

a\ ai | rii , n 2 ) = rii \ ni, n 2 ) ; rii G Z + (132) 

and 

a 2 a 2 | m, n 2 ) = n 2 \ m, n 2 ) ; n 2 e Z + . (133) 



The solution of the recurrence relation which is obtained for the c ni „ 2 's is 

oti a 2 



VW[/(")]I 
with cqo ) fixed as before by normalization, being 



(134) 



oo oo 



\ -1/2 

c oo = ( E E |«i| 2ni |«2| 2n2 ([/(n)]!)- 2 (n 1 !n 2 !)- 1 . (135) 

yni=0 ri2=0 j 

The two-mode f-coherent state can be now defined as 

00 00 a ni a n2 

I oil, oc 2 , f) = c 00 E E / , 1 ■ mi I n i» n2 ) i "i^C; i = l, 2. (136) 
„ 1=0 n 2 =o Vni!n 2 ! [/(n)]! 

It should be remarked that in this form there is a coupling between the two modes, as there is a 
dependence of each of them on the total energy, this interaction between the two modes in general 
is nonlinear. 



Another generalization for two-mode coherent states is of course obtained by means of the 
product of two one-mode f-coherent states, so that there is no interaction between the two modes. 
After defining 

A' t = Oifiim) ; i = 1, 2 (137) 
and finding their eigenstates, we can make the tensor product obtaining 

ai,«2,/i,/2>= ( E E l^in^r" 2 ^!!^!)-^!/!^)]!^)]!^ 21 



i ni=0 ri2=0 



„ 1=0 n 2 =o Vni!n 2 ! [/i(m)]! [/ 2 (n 2 )]! 
In the case f\ = f 2 = 1 , we have the usual two-mode coherent states, namely, 

I a 2 , 1, 1) = e-d^l'^l^ £ E "T== I ^' ™ 2 > • ( 139 ) 

„ 1=0 n 2 =o Vni\n 2 \ 

13 Physical Application of f-coherent States 

We give the examples of how the notions so far discussed might be of some interest in dealing 
with physical problems. With the choice made and repeatedly illustrated, we can continue inter- 
preting | n) as the state containing n quanta, also when it appears in a series representing 
an f-coherent state. Recalling the wide use made in quantum optics of the harmonic oscillator 
formalism, denoting by | n) a state containing n photons, our examples will all be related to 
this interpretation. 

In the usual coherent state | a, 1) , the particle distribution is given by the Poissonian function 

\ n/ \2n 

P 1 , a (n)= l -\e-^ 2 (140) 

with mean photon number 

(n)i = \a\ 2 (141) 

and dispersion 

vi, n =(n 2 ) 1 -(n)j = \a\ 2 . (142) 

Then the ratio of the above quantities cr^ n / (n)i = 1 . 

In the case of f-coherent state the above equations become 

( °° \ a \ 2j \ _ * \ a \ 2n 

p '"<»> = w = (js mJUW) stfwitf (143) 

The mean photon number and dispersion are 

CO 

(n)f = T, nP fM ( 144 ) 

n=0 



and 

oo / oo \ 2 

^/,n = E™ 2p /,»- E™ P /,«H > ( 145 ) 

n=0 \n=0 / 

while their ratio 07, n / (^)/ can be less or greater than 1, producing for a given / either 
sub-Poissonian or super-Poissonian statistics. 

For the case considered in Section 10, where / has zeroes, the photon distribution function 
P On) vanishes for n < n$ and 



N 2 \a 



2n 



P(n) = =M ; n>n . (146) 

n!([/(n)]!) a ' 

The photon distribution in the first of the two generalization for two modes (|128|) is 

|„, |2ni \„, \2no 

P/.c^Cm, na) = c 00 ni!n2 , ([/(n)]!)2 (147) 



and for the second one ( |137|) reads 

P/i,/2,ai,02 = P/i,"i P/2,02 • (148) 

We emphasize the fact that in the latter case the two modes are independent and there is no 
correlation. 

On the contrary, in the previous two-mode generalization a correlation exists between the 
modes: in facts if one gives the interpretation to the f-coherent states as the states related to 
deformation of the annihilation operator due to a specific f-nonlinearity, in the first case, this 
happens through the total energy and then we can now conclude that this particular nonlinearity 
of the field produces a correlation among the modes. In the case of q-oscillators, we obtain for 
the photon distribution in the one-mode q-coherent state 

/ °° | a |2j \ _1 | a |2n 
Pf q ,a( n ) = I E r sinh Aj i| I r sinh An i| ■ (149) 
V=0 i sinh A J - / ^ sinh A J ' 

The property of this distribution is that, for large n (n ^> A -1 ) , the probability to have n 
photons differs essentially from Poissonian distribution due the exponentially decreasing of the 
denominator. The mean photon number (n) q is given for q-nonlinear field by 

sinhAi'n I Z-^i rsinh Ann ' 

(150) 



the second moment by 



r sinh \j 

i sinh A J : / n=0 I sinh A 



\a\ 2 i \ ~ n 2 \a\ 2n 



( n )q — IE r sinh Xj u I E r sinh An l| ' (151) 
\i=0 L sinh X J - / n=0 L sinh A J " 

and the dispersion depends of course on the nonlinearity parameter A . 

The photon distribution demonstrates the fast decreasing of the distribution function in com- 
parison with the Poisson distribution function; the q-nonlinearity makes photon statistics sub- 
Poissonian. The q-nonlinearity and some other types of f-nonlinearities may influence the statis- 
tical properties that can be checked by Hanbury Brown-Twiss-like experiments and the presumed 
detection of sub-Poissonian counting distribution in quantum optics. 



14 Squeezing and Correlation 

Now we will calculate the squeezing and correlation of the quadrature components in the 
introduced f-coherent states. We face, in fact, the problem since the discussed nonlinearity (for 
example, q-nonlinearity) of the field produces the state which is f-coherent state. Then this 
nonlinearity yields the phenomenon of squeezing and correlation of the field (photon) quadrature 
components. It is possible to calculate the dispersion and correlation of the quadratures explicitly. 
To do this we will take advantage of Eq. (|66|), which gives the expression 



f(n + l 
1 



A ; 



/(n + 1) 

Then the quadrature mean values are expressed through 



(a, f | a 1 a, f) = aN^ f 



a 



2d 



n=0 



/(n + l)n!{[/(n)]!} 



I \2 ' 



yelding 



a: 



2» 



a + a* 

(C "' ; ' X ' /} = ~7T Nl > S £ W+ 1) n\ { [f{n)\\ Y 



(a, / I P I a, /) 



a — a 



KsT. 



2n 



^ Q ' ; ^ /(^ + l)^{[/H]!} 2 
The dispersions may be calculated from the relations 

1 2n 



(a, / | a 2 | a, /) = a 2 N 2 aJ £ 



a: 



n==0 



/(n + l)/(n + 2)n!{[/(n)]!} 



| 12 ' 



(a, / | a)a | a, /) = \a\ 2 Nlj ^ 



2(( 



< P{n+l)n\{[f{n)\\y 



as 



1 oo 

(«,/|x 2 |a,/) = -{l + AT2 ;/ [2|«| 2 £ 



2(( 



+ (a 2 + a* 2 )J2 



' P(n+l)n\{[f(n)]\y 
\a\ 2n 



„f f(n+l)f(n + 2)n\{[f(n)]\y 



(152) 



(153) 

(154) 
(155) 

(156) 
(157) 



158) 



1 oo 

(a,f\p 2 \a,f) = -{l + iV^^lal 2 £ 

Z n=0 



2n 



/2(n + l)n!{[/(n)]!p 



Thus, for quadrature dispersion a x = (ft, f \ x 2 \ a, f) — (ft, f \ x \ ft, f) 2 , we have 

a x = ^ + n x a 2 + /i* a* 2 + u x aa*, 



(160) 



where 



*. = *< f E 



In 



< f(n + l)f(n + 2)n\{[f(n)]\y 



N, 



E 



a 



2/i 



< f(n + l)n\{[f(n)W 



(161) 



and 



\a\ 

„t;/ 2 (« + i)«H[/(«)]!}'- 



ft: 



2/i 



^/(n+l)n!{[/(n)]!p 



(162) 



For the other quadrature dispersion a p = (a, f \ p 2 \ ft, f) — (a, f \ p \ a, f) 2 , we have 



a v = - + /i p ft 2 + fi*a* 2 + v p aa* 



(163) 



where 



th = -o N lf\T, 
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f(n + l)f(n + 2)n\{[f(n)]\} 2 
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(164) 



and 



'^ = </ E 
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2ti 



n=0 



f 2 (n + l)n\{[f(n)}\} 2 ^ 



E 



2/t 



' f(n + l)n\{\f(n)]\Y 



(165) 



Depending on the function / (n) the dispersion a x ( cr p ) may become less than 1/2. It means 
squeezing. One can calculate correlation of the quadrature components in f-coherent states as 



a 2 - ft* 2 



0~xp 



N 2 t 

2% aJ 



E 

n=0 



ft 



2/1 



f(n + l)f(n + 2)n\{[f(n)]\y 



— N. 



<*,f 



sr^ \ a \ 
hf(n + l)n\{[f{n)W 



66) 

Then the quadrature correlation coefficient r = a xp / yjo~ x o~ p is not equal to zero. Thus, the 
f-coherent state has the property of being a correlated state [f29]. The invariant a x a p — a 2 p is 
larger than 1/4. Thus, the f-coherent states do not minimize the Schrodinger uncertainty relation 

m, mi. 



15 Deformation of Planck Formula 

We will discuss what physical consequences may be found if the considered f-nonlinearity 
influences the vibrations of the real field mode oscillators like, for example, electromagnetic field 
oscillators or the oscillations of the nuclei in polyatomic molecules. 



First, it will be seen that this nonlinearity changes the specific heat behaviour. To show this, 
we have to find the partition function for a single f-oscillator corresponding to the Hamiltonian 
H = (AA^ + A^A) /2 

n=oo 

Z(T) = £ exp(-(3E n ) , (167) 

?1=0 

where the variable f3 is the inverse temperature T _1 and E n was given in Eq. (^). For the 
evaluation of the quantum partition function for an ensemble of q-oscillators |§, we first note 
that in this case 

E n = — - — - [sinh A(n + 1) + sinh An] ; A = log q , 
smh A 

and obtain that the specific heat decreases for T — > oo as 

Ccc^. (168) 

Thus, the behaviour of the specific heat of the q-oscillator is different from the behaviour of the 
usual oscillator in the high temperature limit. This property may serve as an experimental check 
of the existence of vibrational nonlinearity of the q-oscillator fields. 

The q-deformed Bose distribution can be obtained by the same method and one has HI 



A 2 

(n) = n - (5 — 

D 



\ ((n») - (n)g) + 1 (V) = -n (n%) + (n*) ~ M™ 3 )o 



(169) 



2 V v ,» v , uy 2 
in which no is the usual Bose distribution function and 

(n fc ) = 2 sinh - ^ n fc e " /3(n+1/2) . (170) 

^ n=0 

Calculating the partition function for small q-nonlinearity parameter we have also the following 
q-deformed Planck distribution formula 

1 2 hu e 3h ^ kT + 4e 2n ^ kT + e h ^ kT 

' n ' _ e hu)/kT _ I ~ ~j7p (^ e hu)/kT _ iy ' U'-U 

It means that q-nonlinearity deforms the formula for the mean photon numbers in black body 
radiation p. 

One can write down the high and low temperature approximations for the deformed Planck dis- 
tribution formula For small temperature, the behaviour of the deformed Planck distribution 
differs from the usual one 

(n> - no = -A 2 % . (172) 

For the high temperature, the nonlinear correction to the usual mean photon number also depends 
on temperature 

(n) - n = -A 2 ^ 3 . (173) 

As it was seen, the discussed q-nonlinearity produces a correction to the Planck distribution 
formula (mean oscillator energy) and this may also be subjected to an experimental test. 



As it was suggested in fllQfl , the q-nonlinearity of the field vibrations produces blue shift 
effect which is the effect of the frequency increase with the field intensity. For small nonlinearity 
parameter A and for large number of photons n in a given mode, the relative shift of the light 
frequency is 

5uj X 2 n 2 
~ ~ 2 

This consequence of the possible existence of a q-nonlinearity may be relevant for models of the 
early stage of the Universe. 

Another possible phenomenon related to the q-nonlinearity was considered in |35||, where it 



was shown, that if one deforms the electrostatics equation using the method of deformed creation 
and annihilation operators, a point charge acquires a formfactor due to q-nonlinearity. 



16 Conclusion 

Starting with the example of the harmonic oscillator, we have exhibited a family of associated 
nonlinear systems which are completely integrable, both in classical and quantum physics. 

We have shown that q-nonlinearity, associated with quantum groups, is a subclass of a more 
general class of possible nonlinearity. These aspects, related to the existence of alternative Hamil- 
tonian descriptions for the harmonic oscillator, have been considered with respect to the conse- 
quence for the partition function, again both in the classical and quantum situation. We have 
found that while the partition function for the harmonic oscillator does not depend on the particu- 
lar Hamiltonian, for the nonlinear ones it does depend, giving therefore an experimental possibility 
to select among them. 

A class of states has been considered in the Fock space through the deformation process applied 
to the harmonic oscillator operators. Such states have been described as f-coherent states (or 
nonlinear coherent states), harmonious states and q-coherent states being particular examples of 
them. Their different representation have been constructed, like the Wigner-Moyal and Husimi- 
Kano distributions. It is shown how nonlinear couplings between different modes are easy to 
obtain. 

Keeping unaltered the current physical identification of the Fock states, a possible use is 
presented in the field of quantum optics, obtaining deformed photon distributions and related 
physical quantities. Physical consequences of the deformed vibrations, like the Planck distribution 
deformation, are then reviewed for the q-oscillators, where a blue shift effect exists. The related 



phenomena were studied recently [36-38]. In [36[], the estimation of upper limit of q-nonlinearity 
of the electromagnetic field vibration was done. 

The studied nonlinearities, if they exist, for the electromagnetic field or for the gluons, may 
influence the particle decays, correlations in particle multiplicities, and a change in the Hanbury 
Brown-Twiss experiment results. It would naturally affect the stimulated emission rates and 
hence the radiative equilibrium in the presence of matter. 
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